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Abstract 

The spectral correlation of a chaotic system with spin 1/2 is uni- 
versally described by the GSE (Gaussian Symplectic Ensemble) of 
random matrices in the semiclassical limit. In semiclassical theory, 
the spectral form factor is expressed in terms of the periodic orbits 
and the spin state is simulated by the uniform distribution on a sphere. 
In this paper, instead of the uniform distribution, we introduce Brow- 
nian motion on a sphere to yield the parametric motion of the energy 
levels. As a result, the small time expansion of the form factor is ob- 
tained and found to be in agreement with the prediction of parametric 
random matrices in the transition within the GSE universality class. 
Moreover, by starting the Brownian motion from a point distribution 
on the sphere, we gradually increase the effect of the spin and calcu- 
late the form factor describing the transition from the GOE (Gaussian 
Orthogonal Ensemble) class to the GSE class. 
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1 Introduction 



The universal spectral correlation is one of the most outstanding features of 
quantum systems when the underlying classical dynamics is chaotic [1]. It is 
known that there are universality classes depending on the symmetry of the 
systems. For example, if time reversal invariance is broken, the corresponding 
spectral correlation is reproduced by the GUE (Gaussian Unitary Ensemble) 
of random matrices. On the other hand, the spectral correlation of the 
systems with time reversal invariance depends on the spin. If the system 
is spinless or has an integer spin, the GOE (Gaussian Orthogonal Ensemble) 
gives a precise prediction, while the GSE (Gaussian Symplectic Ensemble) 
applies to a system with a half odd spin. 

In order to explain the universal behaviour from the underlying chaotic 
dynamics, much effort has been paid to establish a semiclassical theory of 
spectral correlations. The spectral form factor K{t) (the Fourier transform 
of the spectral correlation function) is one of the most typical quantities of 
interest. Berry first succeeded in evaluating the leading term in the semi- 
classical r expansion of the spectral form factor [2J. Then Sieber and Richter 
specified the classical orbit pairs which contribute to the second order term[3]. 
More recently Heusler et al. and Miiller et al. extended Sieber and Richter's 
work and calculate the full form of K{t) in agreement with the prediction of 
random matrices [U EJ El El E] • 

In addition to each of the universality classes, the transitions within and 
among them are also of interest. The transitions are described by the spectral 
correlations depending on the transition parameters. It is conjectured that 
such parametric correlations are also reproduced by parametric extensions 
of random matrices[Hl EU]. For the crossover from the GOE class to the 
GUE class, Saito and Nagao invented a scheme to incorporate the transition 
parameters into the semiclassical expansion of -ft'(r)[TT]. Similar schemes can 
also be applied to the transitions within the GUE and GOE classes [13] . 
The agreements with parametric random matrices were in all cases confirmed. 

In this paper, the parametric transition within the GSE symmetry class 
is treated. For that purpose, we shall study the spectral correlation of a 
chaotic system with spin 1/2 by employing the strength of the effective field 
applied to the spin as the parameter. In order to simulate the spin dynamics, 
Brownian motion on the surface of a sphere is introduced. Using semiclassical 
periodic orbit theory, we evaluate the r expansion of the spectral form factor 
up to the third order, so that the agreement with random matrix theory is 
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confirmed. Moreover we study the crossover between a spinless system and 
a system with spin 1/2. We suppose that the Brownian motion starts from a 
point distribution and that a diffusion on the sphere is caused by the increase 
of the coupiing to the effective fieid. As a resuit, the semiciassicai method 
yields the r expansion of the form factor up to the second order. 

The organization of this paper is as follows. In §2, semiciassicai theory of a 
chaotic system with spin 1/2 is developed. Assuming that the spin is coupled 
to a stochastic field, we explain how Brownian motion on a sphere arises. 
Then the leading term in the r expansion of the form factor is evaluated 
by using Berry's diagonal approximation. In §3, a diagrammatic method is 
introduced to calculate the higher order terms in the r expansion. In §4, 
the prediction of random matrix theory is presented and compared with the 
semiciassicai result. In §5, a similar semiciassicai analysis is carried out for 
the crossover from a spinless system to a system with spin 1/2. The last 
section is devoted to a brief summary. 

2 Periodic Orbit Theory for a Chaotic Sys- 
tem with Spin 1/2 

Let us consider the energy level statistics of a bounded quantum system 
with / degrees of freedom. Each phase space point is specified by a vec- 
tor X = (q, p), where / dimensional vectors q and p give the position and 
momentum, respectively. It is assumed that the corresponding classical dy- 
namics is chaotic ( homegeneously hyperbolic and ergo die). Moreover we 
suppose that the system has a spin with a fixed quantum number S. The 
strength of the interaction between the spin and effective field is characterized 
by a parameter 77. 

Let us denote by E the energy of the system. Then, in the semiciassicai 
limit ^ — > 0, the energy level density p{E; Vj) can be written in a decomposed 
form 

p{E-r))^p,,{E)+p^,{E-r)). (2.1) 

Here Pav{E) is the local average of the level density, while Posc{E; r]) gives the 
fiuctuation (oscillation) around the local average. 

The local average of the level density is proportional to the number of 
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Planck cells inside the energy shell: 



Pav(i?) = (25+l)|^, (2.2) 

where the phase space volume with the energy between E and E + AE is 
fl{E)AE. The effective field is assumed to be so weak that Pav(-E') does not 
depend on the parameter t]. 

On the other hand, in order to calculate the fluctuation part posc{E] rj), 
we need to care about the time evolution of the spin. The spin state is 
described by a spinor with 25*+ 1 elements and the spin evolution operator A 
is represented by a (2S'+1) x (25*+!) matrix. We denote such a representation 
matrix evaluated along the periodic orbit 7 by A^(?7). Then, in the leading 
order of the semiclassical approximation, the fluctuation part of the level 
density is written as[6l [HI [15] 



Posc{E;r^) = ^ReY.{tTA,{r^))A,e''-^^^/\ (2.3) 

Here 5*^ is the classical action for the orbital motion, is the stability am- 
plitude (including the Maslov phase) and trA^(r]) is the sum of the diagonal 
elements of A^(?7). 

Now we define the scaled parametric correlation function of the energy 
levels as 



P:.{Ef 

(t^ 

2pav{£;)' 



(2.4) 



Here we introduced averages depicted by the angular brackets (■) over win- 
dows of the center energy E and the scaled energy difference s. The form 
factor, namely the Fourier transform of R{s]r),r)'), is then written as 



/oo 
-00 



I j^^ ^^.rT./n P^^c {e + I;v) Pose (i^-f;^^ 



■ (2.5) 
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Here the angular brackets mean averages over windows of the center energy 
E and the time variable r. Note that r is measured in units of the Heisenberg 
time 

Th = 27rhp,^{E) = {2S + ^^Jfy_, - (2.6) 

It follows from fl2.3p and (12.51) that the form factor is expressed as a double 
sum over periodic orbits 



KiT;v,v') ~ ^ /5:(trA,(r/))(trAy(V))*A^ye^^^^-''^'^/''^ (r 



n \^r-"-'"'~'-' " " " V , 

(2.7) 

where an asterisk stands for a complex conjugate. The periods of the periodic 
orbit 7 and its partner 7' are denoted by Ty and Ty/, respectively. 

In principle, the spin evolution matrix A^(?7) can be calculated from a 
deterministic equation of motion, if the Hamiltonian of the spin is explicitly 
known. However, here we take a simplified strategy based on an assumption 
that the spin evolution parameters undergo Brownian motion on the surface 
of a sphere [T^. The Brownian motion arises when the spin dynamics is 
determined by a stochastic Hamiltonian 

H = r7(h-5), (2.8) 

where 77 is an interaction-strength parameter and S is the spin operator. We 
assume that the components of the effective field 

\V=(h^{t),hy(t),h-A^)) (2.9) 

can be replaced by isotropic Gaussian white noises: denoting the average 
over the noises by the brackets ((■)), we find the correlation 

l,{h,{t)h{f))) = 0, h 

{{hj{t)hj{t'))) = 2D6{t - t') (2.10) 

for j, I = x,y, z. Here isotropy implies that the diffusion constant D does not 
depend on j. 

The time evolution of the spin is described by a (25"+ 1) x {2S+ 1) matrix 
A{t) which satisfies the Schrodinger equation 

in^^A{t) = nA{t), (2.11) 
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where Ti is the matrix representation of the Hamiltonian IK. Note that A(t) 
can be expressed as 

A(t) = exp {i(t){t)SJTi) exp {ie{t)SJh) exp {ii){t)S-Jh) , (2.12) 

where Sx and Sz are (25* + 1) x {2S + 1) matrices representing the x and 
z components of the spin operator S. Thus three Euler angles ip, 9 and 
describe the spin evolution. Let us denote by a segment (with the 

duration T) of the periodic orbit 7. When T coincides with the period, x{T) 
is equated with 7. Along such a segment the spin evolution matrix 

Aj^(r) is evaluated as 

A^(r) = A(T). (2.13) 

Putting fl2.12p into 02.111) . we obtain the Langevin equation for the Euler 
angles 

0/?7 = h^soKpcotO + hyCos(j)cot9 — 
0/rj = — hx cos (j) + hy sin (j), 

ip/rj = — hx sm. (p/ sm 9 — hy cos (f)/ sm. 9. (2-14) 

Then the Fokker-Planck equation 

BP 

— = 7^^DCspP (2.15) 

holds for the p. d.f. (probability distribution function) P{ip, 9, (j)) with the mea- 
sure sin 9dipd9d(j). Here 

^'' = ^9d9''''^d9^^9[W^W'^ """^dm ) ^'-''^ 

is the Laplace-Beltrami operator on the sphere. 

Let us suppose that the Euler angles ip, 9 and are equal to tp',^' and 
(j)', respectively, when the interaction-strength parameter r] is zero. Then the 
solution of the Fokker-Planck equation gives the conditional p. d.f. of the 
Euler angles 

g{^,9,<P;t\^',9',<P') 
= tt t^^Di,^n{i^,9A){DU^\9'A')ye-^^^-''^^'^K 

j=0 m=-j n=-j '^'^'^ 

(2.17) 
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Here -D^^ is Wigner's D function[ 

Di,Jij,9,<p) = e^"^'t'di,je)e'-^, (2.18) 

where 



dLJO) 



(j + m)!(j m)! ^^^m+n^^^^) sin'"-"(g/2)Pjr;;"'"^+"^(cos ( 



^ (j +'^)Kj -^)! 

(2.19) 

with the Jacobi polynomials pI."''^\x). Note that j is an integer or a half odd 
integer (j = 0, 1/2, 1, 3/2, ■ ■ ■ and m, = -j, -j + 1, ■ ■ -J). 

Under the assumption described above, the factor (trA^(77))(trAy (0)) in 
(12.71) with 1]' = can be replaced by the average (( (trA-y(?7))(trAy (0)) )) 
over the Brownian motion. Thus we can write the form factor as 

- ^ (trA,(r/))(trAv(0))* ))A,A;,e^('^-'y^/^6 (r-^^^ 

(2.20) 



We shall evaluate the r expansion of this semiclassical form factor, focusing 
on the systems with spin 5* = 1/2. 

Let us calculate the leading term in the r expansion by using Berry's 
diagonal approximation [2]. In Berry's approximation, one first considers the 
contributions from the pairs of identical periodic orbits (7,7). The spin 
evolution matrix along 7 with S* = 1/2 is given by 

A^(r7) = exp {^1^^ exp {^if^ exp , (2-21) 

where 

are the Pauli matrices. It follows that 

9(1 1 

trA^(r7) = 2cos-cos|-(V' + 0)|. (2.23) 
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The average of the factor (trA^(r7))(trA^(0)) over the Brownian motion can 
be written as 

(( (trA,(r/))(trA,(0)) )) 
= I da;dcu'(trA,(r/))(trA,(O))(7(V^,^,0;T|V'',^',0')Po(^',^',0'), 

(2.24) 

where po is the p.d.f. of the Euler angles at r/ = 0. The integrals are defined 

as 

dcu = / dij d6 d0sin^, 
JO Jo Jo 

du' = / dtp' de' / d0'sin^' (2.25) 

JO JO JO 



and T = Ty is the period of 7. 

For the transition within the GSE universality class (the GSE to GSE 
transition), we employ the uniform "initial distribution" 

Po{i^,OA) = ^, (2.26) 

since it yields the spectral form factor of the GSE class [H [HI [151 [IS]- The 
uniform distribution at 77 = implies that the spin is under the influence of 
additional interactions apart from the interaction described by (12. 8p . Putting 
([22S]) into dOl, we obtain 

(( (trA,(r/))(trA,(0)) )) 
= ^ / du;du;'{tTA,{r])){tTA,{0))g{ij, 9, T|^', 9', </>') 

= ^ / duduj' {/^;;2,i/2(^, 0, 0) + D%^,,^,i^, 9, (/>)}* 

(2.27) 



Therefore, using the definition (I2.17P of g and the orthogonality relation 

327r^ 
2j + l 



D^„,„(^, ^, 0)} „,(V^, ^, 0) = ——6,,,Sm,m'Sn,n', (2.28) 
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we can readily find 

(((trA,(r^))(trA,(0))))=e"(^/^)"^ (2.29) 

with 

a = rfD. (2.30) 

Here the interaction-strength parameter rj is scaled so that aT remains finite 
in the semiclassical limit /i — > 0. In order to take a step further, we need 
Hannay and Ozorio de Almeida (HOdA)'s sum rule[ 



which results from the ergodicity of the system. Using this sum rule, we find 
the contribution to the form factor as 

i^(7,7)(^;^,0) = ^(|ElAl'^(^-^))(((trA,(r/))(trA,(0)))) 

= re-(=^/^)'^^. (2.32) 

The second contribution to Berry's diagonal approximation comes from the 
pairs (7,7), where a bar denotes time reversal. Noting 



A,(r^) = {A,(r^)} 



-1 



= [li (2.33) 

where {A^(r7)}'^ is the transpose of A-y(?7), we find 

trA;^(r7) = trA^(r7). (2.34) 
Therefore we can similarly obtain a contribution 

^(,,^)(r;r/,0) = A.(|5]|A,|^5(r-^)^(((trA,(r/))(trA,(0)))) 

= re-(=^/^)'^^. (2.35) 
Thus the total sum of the contributions to the diagonal approximation is 

i^diag(r) = K(^,,,){t; V, 0) + K^,,^){t; v, 0) = 2re-(=^/^)'^^. (2.36) 
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3 OfF-diagonal Contributions 



We are now in a position to calculate the off-diagonal contributions, restrict- 
ing ourselves to the systems with two degrees of freedom (/ = 2). Encounters 
of periodic orbits play the major role in identifying the leading terms. An 
encounter is a set of orbit segments which come close to each other in the 
phase space. Long periodic orbits have encounters of the order of the Ehren- 
fest time Te- In the semiclassical limit h —>■ 0, Te logarithmically diverges. 
However, as the period T is of the order of the Heisenberg time Th, which 
more rapidly diverges, Te remains vanishingly small compared with the pe- 
riod. Therefore the periodic orbit mostly goes along loops in the phase space 
and occasionally visit encounters. As the leading terms are expected to result 
from the periodic orbit pairs (7, 7') which are close to each other or mutually 
almost time reversed, we can suppose that 7' is almost identical to 7 or 7 on 
the loops but differently connected in the encounters. 

Let us consider such a periodic orbit pair a = (7, 7') in the phase space. 
Within each encounter, a Poincare section V orthogonal to the orbit 7 can 
be introduced. Suppose that 7 pierces V within the r-th encounter. If Ir 
segments of 7 are contained in the r-th encounter, there are piercing points 
on V. The displacement 5x between such points can be spanned as 5x = 
sis + uiu- Here pairwise normalised vectors ig and iu have directions along 
the stable and unstable manifolds, respectively. Therefore, if one reference 
piercing point is chosen as the origin, each of other piercing points is identified 
by a coordinate pair {s,u). As a result, if 7 has L loops and V encounters, 
Y.r=ii^r — 1) = L — V coordinate pairs {sj,Uj), j = 1,2, ■ ■ ■ , L — V are 
necessary to identify the piercing points of 7. 

Let us denote by Tj the duration on the j-th loop and by the duration 
of the r-th encounter. Then the total duration of the encounters is 

V 

ta = ^rtr (3-1) 
r=l 

and the period is 

T = X:T,+t„. (3.2) 
Ergodicity can be employed to estimate the number of encounters as[5l [6|, [TJ 

m 

/rT — ta rT — ta—Tl rT — ta—Tl—T2 71,-2 

duds/ dTi / dT2---/ dTi_ig„ (3.3) 

JO Jo Jo 



10 



where the integration measures are given by 



L-V L-V 

du = n dwj, ds = n dsj (3.4) 

and 

T 

Here is a combinatorial factor chosen such that overcountings are avoided. 

Now the contribution to the form factor from the orbit pair a = (7, 7') 
and its counterpart (7,7') can be readily derived. Referring to fl2.20l) and 
(13.31) and taking account of (12.341) . we find that such a contribution is 

/rT—ta rT-ta-Ti rT-ta-Ti—T2 

duds / dTi / dT2 ■ ■ ■ / dTi_i 

Jo JO Jo 

X Q„(( (trA^(r7))(trAy(0)) ))e^^^/^ (3.6) 
where the action difference AS* = — Sy is given by AS* = J2j=i %Sj[5l, EJ 

In order to obtain a semiclassical result, we need to expand Ka{T) in t^'s 
and extract the terms in which all t^s mutually cancel. Since extra factors 
h appear or rapid oscillations take place in the limit ^ ^ 0, the other terms 
should be neglected [5l[6l [7]. The off-diagonal contribution to the semiclassical 
form factor is thus derived as 



tl=t2 = ---=tv=0 

(3.7) 



where 



$(ti,t2,---,ty) 

T — ta pT — ta—Ti rT — ta— T1-T2 Tl-2 

dTi / dT2 ■ ■ ■ / dTL-i 



X (((trA,(r^))VAy(0)))). 

(3.8) 



3.1 Sieber-Richter Term 

In this and next subsection we consider the r expansion of the above formula 
(13.71) . Mathematica was used to assist the computations. Each term of (13.71) 
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is of order r" with n = L — V + 1. Let us first consider the second order 
term {n = 2). The relevant pairs a = (7,7') has two loops {L = 2) and 
one encounter {V = 1). Such periodic orbit pairs were identified by Sieber 
and Richter and thus called SR (Sieber- Richter) pairs[3]. An SR pair is 
schematically depicted in Figure 1. 




Figure 1: The periodic orbit pair contributing to the second order term 

In Figure 1, Li and L2 are loops and Ei is an encounter. In the encounter, 
7 and 7' are depicted by solid curves and dashed lines, respectively, and each 
arrow shows the direction of the motion. We can symbolically write the 
periodic orbits as 

7 = ^iLs^iLi, 7' = E[L'2E[L[, (3.9) 
so that the spin evolution matrices are 

= (AeX'^l^Ae^Al,, Ay = {Ae'X\^lO~'^e[Al',. (3.10) 

A spin evolution matrix A^ along a segment x '^^ ^ periodic orbit is given by 
(I2.13P and can be expressed as 



A^ = exp [(I>x2^zj exp (0x2^=^) exp (i^x2^zj (3.11) 

in terms of a set of the Euler angles = (V^^, 6'^, 0^). The Pauli matrices 
a^c and are defined in (12.221) . The corresponding integral over the Euler 
angles is defined as 

//'47r rn rin 
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Moreover we denote the durations of Li, L2 and Ei by Ti, T2 and ti, respec- 
tively. Using the above notations, we evaluate the average of (trA^(r7))(trA^/(0)) 
as 



(327 



(( (trA,(r^))(trAy(0)) )) 

J dUL^UL^duEi j dcu^/dt^L/^dc^ij/ 

X tr((AsJ-iAi,As,AiJtr((A^j)-i(A^,)-iAsjA^;) 
X g{ujL^;Ti\ujL,^)g{ujL2]T -Ti - 2ti\ujL'^)g{ujE^\ti\ujE[) 

^ lg-{3/4)aT ^^m)aU _ 3e-{l/2)aii j ^ (3^ ^3) 

Due to the equivalence of the segments Ei and ^i, we need to choose the 
combinatorial factor Na as N^^ = 2^]. Consequently we find the contribu- 
tion from the SR pairs to the form factor as 



so that 



i^SR r = -— — $ ti 

A'^SR Oti 



2rV(3/4)«r(l + ^ar). (3.15) 



3.2 Third Order Term 



Next we consider the third order term (ra = L — ^ + 1 = 3). It is known that 
the periodic orbit pairs contributing to the third order term are classified 
into five types: aas, api, ppi, ac and pc[4j. These five types are depicted in 
Figure 2. 

As is seen from Figure 2, each of aas, api and ppi orbit pairs has four 
loops (L = 4) and two encounters {V = 2). The durations of the loops 
Lj (j = 1,2,3,4) and the encounters Ei (/ = 1,2) are denoted by Tj and 
ti, respectively. The combinatorial factors are known to be given by 
iVaa. = 2, A^api = 2 and A^ppi = 40]. 

On the other hand, each of ac and pc orbit pairs has three loops {L = 
3) and one encounter {V = 1). The times elapsed on the loops Lj (j = 
1, 2, 3) and on the encounter Ei are denoted by Tj and ti, respectively. The 
combinatorial factors Na are A"ac = 1 and A^pc = 3[4J. 

In the following we calculate the contribution to the form factor Ka{T) 
from each of the five types: a = aas, api, ppi, ac and pc. 
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aas 




ac pc 

Figure 2: The periodic orbit pairs contributing to the third order term 

(1) aas orbit pairs (A^aas = 2) 

(( (trA,(ry))(trAy(0)) )) 
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16 

Therefore 



(327r2)6 

X / dujL'dujL'dujL'dujL'dujE'dujE' 

I -^1 -^2 -^3 ^4 ^1 -'^2 

X tr{AEAL2^E2^Ls{^E2)~^^u{^E,)~^^L,) 

X tr(A^,(Ai,)-^A^,A^,(A^,)-^(Ai,)-^(A^,)-^Ai,) 

X g{ujLi]Ti\ujL'^)g{ujL2]T2\ujL'^)g{ujL^] Tz\ujl'^) 

X ^K,; T-2ti- 2t2 -Ti-T2- Ts\ujlO 

X 5r(a;£;i;ii|u;E/)5r(a;E2;i2|<^£;^) 

J_g-(3/4)aT J^g(3/2)ati _ 3g-(l/2)ati^ ^g{3/2)at2 _ 3g-(l/2)at2^ _ (^g jg^ 



$(^1,^2) 

{T-2ti - 2t2)-'^ ^_(3/4)„-r j^g(3/2)ati _ 3g-(l/2)ati j j^g(3/2)at2 _ 3g-(l/2)at2j ^ 

(3.17) 



SO that 



ti=i2=0 

= 4T3e-(3/4)«^ \l + \aT+ ^{aTf] ■ (3.18) 
(2) api orbit pairs (A^api — 2) 



(327r2)6 

X 



(( (trA,(ry))(trAy(0)) )) 

j dUL-^dUL^dUL^dUL^dUEidUE^ 

j dujL'^dujL'^dujL'^dujL'd^E[du:E'^ 

X tr ( A^^ A^^, As, Ai, ( As J -1 Ai J 

X tr(AE(Ai/AE^(As^)-i(AE()-'(As'J-'AE^Ai.) 

X ^ (c^L, ; Ti I c^L'J ^ (c^L2 ; ^2 1 ) ^ ( ; ^3 1 Wl^, ) 
X ^(^i,; T - 2^1 - 2^2 -Ti-T2- Ts\iUL') 
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16 

It follows that 



X g{u;E^;ti\u;E'Jg{ujE2;t2\uJE:^) 



i^api(r) = 4r^e-(^/^)'^^ {l + ^aT + ^{aTf] . (3.20) 

(3) ppi orbit pairs (iVppi = 4) 

(( (trA,(^))(trAy(0)) )) 
= (^^2^2y J dujLidu;L2dujL3dujLAujEidujE2 

X / 

diOL'j^ duji,'^ dujij'^ duj^'^ dc*.^^^ dujEi^ 

X tr{AE,ALAEALs^EALAEAL,) 

X tr(AB.Ai.AE^Ai.AEiAi.AE.Ai.) 

X g(u;L,;Ti\u;L[)g{u;L2]T2\ujL'Jg{u;Ls; TsIljl'^) 

X g{uju; T-2h- 2t2 -T1-T2- T^\oJlO 

X g{uEi]ti\uE[)g{<^E2]h\^E'^) 

^ J_g-{3/4)aT J^g(3/2)ati _ 3g-(l/2)ati j |^p{3/2)at2 _ 3g-(l/2)at2^ _ (3 21) 

It follows that 

Xpp^r) = 2r^e-(^/^)«^ {l + + |(«^)'} ■ (3-22) 

(4) ac orbit pairs (A^ac = 1) 

(( (trA,(r;))(trAy(0)) )) 

= ^3^^2)4 / dUL^dUL^dULzdUEi j dUL'dUL'^dUL'^dUE[ 

X tr(As,Ai,(A^J-iAz.,A^,Ai3)tr(As.(A^,)-i(Asj)-i(A^,)-iAEjA^.) 
X g{ujL^;Ti\ujL'^)g{ujL2]T2\uJL'^)g{'^L3]T - 3ti - Ti - T2\uJl'^) 

X g{uJE^,ti\uJE'^) 
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Therefore 

8 ^ ^ 

so that 



(5) pc orbit pairs (A^pc = 3) 



ti=0 



4 

It follows that 



(( (trA,(77))(trAy(0)) )) 
= ^g^^2)4 / dc^LidcuLadcuLsdcUiJi ^ dcjz.; dw^^ dt^L^ du;£;j 

X tr ( Ai^, Al, Ae, Al, Ae, Ai3 ) tr ( A^j. A^, A^^j A^, A^^j A^, ) 
X g{.^L^]Ti\uL:^)g{ujL^]T2\ujL09{.^L-i,T - ?>ti -Ti -Tslc^L^) 
X 5((u;£;i;ti|u;£;j) 

ie-(3/4)'^T |^2e-(3/2)-*i _ ^m)at^^^ . (3.26) 



Kpe(r) = -2T^e-'^'/^)-^ [^ + \(^T]. (3.27) 



Putting the above results together, we obtain the third order contribution to 
the form factor 

i^ard (r) = i^aas (r ) + i^api (t ) + iTppi (t) + Kac (t) + Kpe (t ) 

= 2rV(=^/^)'^^|l + ^aT + ^(aT)^|. (3.28) 

Hence the semiclassical form factor up to the third order is calculated from 
(EJSD, (13151) and (13:281) as 

irsc(T) = i^diag(r) + irsR(r) + ir3rf(r) 

3 _\ r. 3 _ 15, 



2^g-(3/4)aT 



l+{l + \aT)r+[l + \aT + f^{aTfy 



(3.29) 
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4 Parametric Random Matrix Theory 



Parametric random matrix theory was originally invented by Dyson [20]. The 
quantum Hamiltonian of a time reversal invariant system with spin 1/2 is 
simulated by an x self-dual real quaternion random matrix H. It is 
assumed to be a sum of a self-dual real quaternion matrix Hq and a Gaussian 
random perturbation: the p.d.f. of H is given by 



P{H-a\Ho) dH ocexp 



Tr{(if-e-/7o)'} 



1 -e- 



-2(T 



dH 



with 



TV N 3 

dH=YldH,,Y[Y[dH]t 

j=l j<l k=0 



{k) 



(4.1) 



(4.2) 



Here ifjf'* is the k-th component of the real quaternion Hji. We are interested 
in the parametric motion of the matrix H depending on the fictitious time 
parameter a. 

Let us write the eigenvalues of the self-dual real quaternion matrices H 



and Hq as xi, X2,---,xn and yi, 1/2, ■ 
Fokker-Planck equation 

dp _^ d 



■,yN, respectively. Dyson derived the 



dW 1 dp 

OXj 



4 dxj 



with 



N 



w 



N 

^log|Xj 



(4.3) 



(4.4) 



for the p.d.f. p of the eigenvalues of H. 
We denote by 



G{xi,X2, 



.XN;cr\yi,y2, 



• Vn) 



(4.5) 



the Green function solution of the Fokker-Planck equation (14. 3p . Namely, G 



with the measure D^i dxj gives the p.d.f. of the eigenvalues of if at cr under 



the condition that Xj = yj {j 



the Green function is given by the p.d.f. of the GSE eigenvalues 
G{xi,X2, ■■■,xn; oo|?/i, ?/2, ■ ■ ■ , Vn) = PGSE(a:i, 2:2, ■ ■ • , Xn), 



1, 2, • ■ ■ , N) at (T = 0. The limit a — > 00 of 

(4.6) 
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where 

Pgse{xi,X2,- ■ ■ oce"^^ . (4.7) 

Let us choose the initial matrix Hq as a GSE random matrix. Then 
the transition within the GSE symmetry class (the GSE to GSE transition) 
is realized. We define the dynamical (density-density) correlation function 
describing the correlation between the eigenvalues of H and Hq as 

K{x-a\y)=N^^-^-^, (4.8) 

-'0 

where 

dx2 / dx3 ■ ■ ■ / dxN / d?/2 / d?/3 ■ ■ ■ / dyN 

-OO J —CO J —oo J —oo J —oo J —oo 

X G{xi,X2, ■■■,Xn] cr\yi, 2/2, ■ ■ ■ , 2/iv)PGSE(l/l, 2/2, " " " , 2/iv) 

(4.9) 

and 

Jo= / dx / d2/J(x;cr|2/). (4.10) 

J —oo J —oo 

The asymptotic limit — * oo of the dynamical correlation function was 
evaluated by the method of supersymmetry [21j . It can also be derived by 
using the properties of the Jack symmetric polynomials [22]. Let us note 
that the asymptotic eigenvalue density at ^/2Nz (— 1 < z < 1) is given by 
p = \J2N{1 — z^)l 71. In terms of the new scaled variables c, X and Y defined 
as 

a = c/(7rV), x = V2Nz + {X/p), y = V2Nz + {Y/ p), (4.11) 

one obtains the asymptotic limit 
k(x: cr\y) 



/oo /■! /•! 

du J dvi J dv2 



[v^ -1){U- V1V2) 



2 



2 A J-i 'J-i ' {2uviV2 ~ - {viy - {V2y + IV 
^ ^-cW+iv,f+{v2?-2MHv2f-i}cos{2n^{u-viV2)} (4.12) 

with ^ = X — Y. The Fourier transform of the asymptotic limit 

/oo 
dee*'""«p(e;c) (4.13) 
-00 
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gives the definition of the form factor. It can be written as 



^rm(t) 



r-2 .1 



Jl-T J(l- 



dvo 



{1-t)/vi 



X 



{V1V2 + t)2 - 1 



{2{VIV2 + t)viV2 - {ViV2 + rf - - {V2f + 1}2 



^ ^-c{(viV2+Tf+(viy+{v2f-2{vif{v2f-l} (4-14) 

for < T < 1. In order to derive the r expansion of KnuiT), we introduce 
new integration variables Si and S2 by 



T T 
Xl — 1 — —Si, A1A2 = 1 — —S2- 



(4.15) 



Then we find 



Krm{t) = ^ ^ dsi ^ ds2 exp |-2A (l - ^^2) - Arj 



X exp < 



At 



X 



2 ' 



(4.16) 

where A = cr. Thus we can readily calculate the r expansion (with fixed A) 
from the Taylor expansion of the integrand as 

KMr) = {4 + (2 + 4A)r + (l + 2A + yA^) + ■ ■ j . (4.17) 

In order to compare this result with the semiclassical formula, we need 
to take account of the Kramers degeneracy, which means that all the eigen- 
values have multiplicity two due to time reversal symmetry. Inclusion of the 
degeneracy yields a modified form factor 



Knuir) = 2XRM(2r) 

= 2re-^^ {1 + (1 + 2A)r + (l + 2A + ^A^) + • • j . 



(4.18) 
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This is in agreement with the semiclassical formula fl3.29l) up to the third 
order with an identification A = (3/8)aT. 



5 The GOE to GSE Transition 

If the spin evolution operator is represented by an identity matrix, the system 
is effectively spinless and the resulting spectral correlation belongs to the 
GOE universality class. Therefore, the crossover from the GOE class to the 
GSE class can be treated by introducing 

PolV-, 0) = 5{iP)5{cos e - 1)5(0) (5.1) 

as the "initial distribution" instead of fl2.26p . In this section we investigate 
the GOE to GSE transition, focusing on the form factor i^(r, 77,77), where rj' 
is equated with 77. 

As before, due to the relation fl2.34p . the contributions from the pairs 
(7, 7') and (7, 7') are equal. Therefore, in order to calculate the form factor in 
the diagonal approximation, it suffices to treat the pairs (7,7). The average 
over the Brownian motion can be evaluated as 

(( (trA,(7/))2 )) 
= j dujduj'{tT/\,{r,) fg{,p, e, 0; T|V^', 9', 0')Po(V^', 9', 0'), 

= J duj{tTA,i7^)fg{^, e, 0; r|0, 0, 0). (5.2) 

Noting 

(trA,(77))2 = D%i^P, 9, 0) + D\_,iij, e, 0) + Dl.ii;, 9, 0) + Dl.O;, 9, 0) 

(5.3) 

and the orthogonality relation (12.281) . we can readily find 

(((trA,(7^))2)) = l + 3e-2'^^. (5.4) 

Then, using the HOdA sum rule fl2.3ip . we find the contribution to the form 
factor 

Ki,,,){T;v,v) = ^(ElAl''^(^-^))(((trA,(7?))^)) 

= r(l + 3e-2'^^), (5.5) 
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so that the diagonal term arising from the pairs (7, 7) and (7, 7) is 

i^diag(r) = i^(^,-,)(r; r/, r/) + K(^,^)(r; r], r]) = 2r(l + Se-^"^). (5.6) 

Let us next consider the second order term. As before, it can be evaluated 
from the Sieber-Richter pair (7,7') in Figure 1. We compute the average of 
(trA^(?7))(trAy (77)) over the Brownian motion as 

(( (trA,(r/))(trAy(r/)) )) 

= J duJLiduJL2(i^Ei 

X g{iJL,; Ti|0, 0, 0)g{u;L,; T - Ti - 2ti|0, 0, 0)g{u;E,; h\0, 0, 0) 

^ — 1 + ^„-2aT+4ati ^ -2aT+2aTi+4ati , ^ _2aTi ^\ 

2 2 2 2 ■ ^ ' ' 

Then we can evaluate the contribution to the form factor 

= 2r2{l + (6aT-9)e-2'^^}. (5.8) 

Thus we obtain the semiclassical form factor up to the second order 

Kscir) = i^diag(r)+KsR(r) 

= 2r(l + 3e-2'^^) + 2r2{l + (6aT-9)e-2'^^}. (5.9) 

A random matrix model of the GOE to GSE transition was already for- 
mulated in 1^ . However, as far as the authors know, an asymptotic 
formula to be compared with the above result (15.91) has not been worked out. 
Therefore it can be regarded as a conjecture for one of the open problems in 
random matrix theory. 

The corresponding random matrix model can be formulated by using 
Dyson's p.d.f. (14.11) . Here we need to suppose that the initial matrix Hq is a 
GOE random matrix. Namely, the matrix elements of Hq only have the 0-th 
components and the p.d.f. of Hq is 

PGOE{Ho)dHo oc e-(^/2)Tt(//o)^di/o (5.10) 
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TV TV 

dHo = f[d{Ho)j,f[d{Ho)ji. (5.11) 

It is well known that the form factor of the GOE eigenvalues is expanded as 

KooEir) = 2t - 2t^ + ■ ■ ■ . (5.12) 

Considering the Kramers degeneracy, one modifies it into 

KcoEir) = 2irGOE(2r) = 8r - IQr^ + ■ ■ ■ , (5.13) 

which is in agreement with the corresponding case a = of the semiclassical 
result ra . 

6 Summary 

In this paper, the parametric spectral correlation of a chaotic system with 
spin 1/2 was studied. The parameter was chosen to be the strength of the 
effective field applied to the spin. Using the semiclassical periodic orbit 
theory for the orbital motion and simulating the spin dynamics by Brownian 
motion on a sphere, we evaluated the parameter-dependent spectral form 
factor Ksc{t)- The r expansion of ii'sc(T) was found to be in agreement 
with the prediction of random matrix theory up to the third order. Moreover 
a crossover from a spinless system was investigated and the r expansion of 
the corresponding form factor was calculated up to the second order. 
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